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ABSTRACT 

We present an analysis of the kinematics and excitation of the warm ionized gas in 
two obscured, powerful quasars at z^3.5 from the SWIRE survey, SWIRE J022513.90- 
043419.9 and SWIRE J022550.67-042142, based on imaging spectroscopy on the VLT. 
Line ratios in both targets are consistent with luminous narrow-line regions of AGN. 
SWIRE J022550.67-042142 has very broad (FWHM=2000 km s'^), spatially compact 
[OIII] hue emission. SWIRE J022513.90-043419.9 is spatially resolved, has complex 
line profiles of II/3 and [OIII], including broad wings with blueshifts of up to —1500 km 
relative to the narrow [OIII]A5007 component, and widths of up to FWIIM=5000 
km s~^. Estimating the systemic redshift from the narrow H/3 line, as is standard for 
AGN host galaxies, implies that a significant fraction of the molecular gas is blueshifted 
relative to the systemic velocity. Thus the molecular gas could be participating in the 
outflow. Significant fractions of the ionized and molecular gas reach velocities greater 
than the escape velocity. We compare empirical and modeling constraints for different 
energy injection mechanisms, such as merging, star formation, and momentum-driven 
AGN winds. We argue that the radio source is the most likely culprit, in spite of 
the sources rather modest radio power of 10^^ W Hz^^. Such a radio power is not 
uncommon for intense starburst galaxies at z^2. We disucss these results in light of 
the co-evolution of AGN and their host galaxy. 

Key words: galaxies: formation, galaxies: high-redshift, quasars: emission lines, 
galaxies: kinematics and dynamics 



1 INTRODUCTION 

Massive galaxies provide important constraints for our un- 
derstanding of galaxy evolution. Observational and theoret- 
ical progress has led to a reasonably clear picture: Massive 
galaxies formed most of their stars in relatively short, in- 
tense bursts at high redshift, that were triggered by ma- 
jor, gas-rich mergers, and regulated by the energy injection 
from star formation and active galactic nuclei (e.g., Hop- 
kins et al.||2006[ [Granato et al.1|2006| [Chakrabarti et al. 
Narayanan et al.||2008 |. This basic scenario is simple 



2008 



and elegant, and broadly consistent with ensemble studies 
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of statistical samples of galaxies with intense star forma- 
tion and powerful AGN. However, given the complexity of 
these processes, neither ensemble studies nor models alone 
can show conclusively if the different components of this 
model - merger, star formation and AGN activity - interact 
as postulated, and if this has the predicted impact on the 
interstellar medium. This is particularly true for the role of 
the energy injection by the AGN, for which we do not yet 
have a good physical understanding. 

The most direct way of overcoming these limitations 
is through detailed observations of the gas kinematics and 
energetics in massive galaxies that are in key phases along 
this sequence. Ideally, such an approach must focus on an tn 
situ study of massive galaxies during their major phase of 
growth - hence we must observe galaxies at high redshift. 

Here we present an analysis of the warm ionized in- 
terstellar medium in two such galaxies, which is based on 
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deep rest-frame optical integral-field spectroscopy obtained 
with the Very Large Telescope of ESO. SWIRE J022513.90- 
043419.9 and SWIRE J022550. 67-042142. 2 (SW022513 and 
SW022550 hereafter) are two obscured quasars at z^3.5 
( Polletta et al.pOOS ) and are the two most luminous 24^m 



emitters at high redshift in the SWIRE survey (Lonsdale 
et al.|2003 l. We also include the results of a recent CO (4-3) 



analysis into our discussion, which has been presented by 



Polletta et al. (2011 



Polletta et al. ( |2008[ ) presented a detailed analysis of 



the multi-wavelength photometric properties of these galax- 
ies. Both galaxies host powerful starbursts {£, ~ W^^Cq) 
and luminous, obscured quasars with bolometric luminosi- 
ties > 10^"^ Cq. Centimeter radio observations suggest the 
presence of moderately powerful radio sources of order 10^^ 
W Hz~^ at 1.4 GHz in both targets, which appear too pow- 
erful to be entirely powered by star formation. The rest-UV 
spectrum of SW022550 shows bright, high-ionization emis- 
sion lines, in particular NVA1240, and weak UV continuum 
emission, which are the typical signatures of luminous, ob- 
scured quasars ( [Polletta et"aL]|2008|). SW 022513 is well de- 
tected in the X-ray (Pierre et al. 20071 with a luminosity, 



£(2-10 keV) = 8.2 x 10*" erg s"^ in the hard X-ray that 
is 7 times brigher than the soft X-ray luminosity at 0.5-2 
keV. It is thus a luminous absorbed X-ray source, e.g. a fac- 
tor 3 brighter than the archetypal type-2 quasar CDFS-202 
at z=3.700 ( [Norman et al.|20 02). SW022513 is most proba- 
bly Compton-thick since its hard X-ray band luminosity is 
about only 1/50 of that measured around 24/im. 



Polletta et al. ( 2008 1 also present low-resolution ISAAC 



longslit spectroscopy in the near-infrared for both targets, 
and rest-frame UV spectroscopy for SW022550, showing 
that both sources are luminous line emitters. Recent IRAM 
Plateau de Bure millimeter interferometry of the CO (4-3) 
line shows that both galaxies are luminous CO line emitters 
with Ceo ~ 5 X 1O^*^£0 K km s~^ pc^ ea ch, correspond- 
ing to 4 X 10^" M0 of cold molecular gas (Polletta et al 



20111. SW022550 has a broad, double-horned CO line pro- 



file, whereas the CO(4-3) line in SW022513 is also broad, 
FWHM=1000 km s~^, and featureless, and is not very well 
fit with a Gaussian profile. Short gas consumption timescales 
suggest the galaxies may be near the end of their epoch 
of intense star formation, which is when galaxy evolution 
models postulate the impact of the AGN should be greatest 
( Springel et al.„2005l [Narayanan et al.|2008 l. 

Targeting particularly powerful quasars is important to 
isolate the impact of the AGN from that of the starburst 
(or potentially a merger). By contrasting our targets and 
galaxies that have equally intense star formation, but less 
powerful AGN we may hope to identify the unique signa- 
tures of the AGN. In addition we may set upper limits on 
what impact a luminous, obscured AGN may possibly have 
on the gas kinematics of its host galaxy. Both aspects are 
important to quantify the impact of AGN on the evolution 
of their host galaxies, as they may serve as benchmarks to 
develop analyses of less powerful, more frequent forms of 
feedback later on - and to infer whether these are indeed 
observationally feasible with present-day instruments. 

Both galaxies have kpc-sized narrow-line regions 
(NLRs) that have strongly disturbed gas kinematics and for- 
bidden line emission with widths of up to few 1000 km s~^. 
SW022513 has complex line profiles with broad, blueshifted 



components seen in [OIII]AA4959,5007 and H/3. These are 
qualitatively similar, but broader than the blue wing previ- 
ously found by Alexander et al. (2010) in the submillimeter- 
selected quasar SMMJ1237-f 6203. Our analysis focuses on 
the gas kinematics, including the ionized and molecular gas 



studied by Polletta et al. (20111. We compare with the en- 



ergy output of SW022513 and SW022550 through AGN and 
star formation and discuss what physical processes may con- 
tribute to driving such an outflow. 

Throughout our analysis we adopt a flat cosmology 
where Ho = 70 km s"^ Mpc"\ Qa^O.!, Qm = 0.3. With this 
cosmology, the luminosity distance to SW022513 at z=3.42 
is 29.7 Gpc. 1 arcsec corresponds to a projected distance 
of 7.4 kpc. SW022550 is at a luminosity distance Dl =34.4 
Gpc, and 1 arcsec corresponds to a projected distance of 7.0 
kpc. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed both targets with the SINFONI integral-fleld 
spectrograph on the VLT. SINFONI is an image sheer which 
gives a held of view of 8" x 8" with a pixel scale of 250 mas 
in the seeing-limited mode. We used the H-l-K grating which 
covers the near-infrared H and K bands simultaneously from 
1.45/^m to 2.45/im at a spectral resolving power of R~1500 
(~200 km s"i). 

Both targets were observed in October and November 
2009 under good and stable atmospheric conditions. The 
program was carried out in service mode under Program ID 
384.B-0161. We obtained 3 hrs of on-source observing time 
per target, split into individual exposures of 300 seconds to 
allow for a good sky subtraction. Since we expected that 
both sources were smaller than the field of view we adopted 
a dither pattern where the sources fall onto different parts of 
the detector, but within the field of view at all times during 
observations. This allows to use one exposure to subtract 
the sky from the subsequent exposure, and makes taking 
dedicated sky frames unnessesary. 

Data reduction relies on the standard IRAF tools to 



reduce longslit spectra (Tody 19931, which we modified to 



meet the special requirements of integral-field spectroscopy. 
Our data reduction has been extensively described elsewhere 
(e.g., [Nesvadba et al.|[2006[ [2007| 2008). The absolute flux 



scale was determined from standard star observations at the 
end of each one-hour long observing sequence, and we used 
the same stars to measure the seeing, which typically has a 
fuU-width-at-half-maximum of 0.5" to 0.6". 



3 RESULTS 

3.1 Continuum emission 

We detect continuum emission at the position expected from 
the IRAC 3.6fim imaging in both galaxies. The spectral 
coverage of our H-j-K data corresponds to 3300A to 5500A 
and to 2900A to 5000A in the rest-frame of SW022513 and 
SW022550, respectively. The continuum is relatively faint 
in both galaxies, but nonetheless detected at a significance 
of 1— 2cr per pixel in each spectral bin of SW022513 and 
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SW022550. To obtain a more robust measurement, we col- 
lapsed both cubes at wavelengths without strong line emis- 
sion, and detected the continuum at significances of >3cr and 
up to 14(T per spatial pixel in the collapsed image over the 
area of one PSF. We did not attempt to constrain the over- 
all spectral shape of the continuum, however notice that it 
seems to have a blue slope in both targets broadly consistent 
with originating from the AGN. In this case, the continuum 
could be either produced by scattered or direct light from 
the AGN, or represent nebular continuum (e.g., Vernet et al. 



2001 1, A more detailed discussion of the continuum emission 



is beyond the scope of this paper. 

The continuum of SW022550 is compact, in SW022513 
it is marginally extended along the North-South axis at ~ Scr 
significance (per spatial pixel, see Figure [2|. We also detect 
very faint continuum emission associated with the faint line 
emission north from the nucleus in SW022513. 



3.2 Emission-line gas 

3.2.1 SW022513 

We show the integrated spectrum of SW022513 in Figure [l] 
Line properties are broadly consistent with those found by 



PoUetta et al. ( 2008 1 with ISAAC longslit spectroscopy at 



3x lower spectral resolution and lower signal-to- noise ratios. 
The lines are spectrally well resolved with typical widths 
of FWHM^700 km s"\ Due to the broad width of the 
lines we did not resolve the individual components of the 
[OII]AA3726,3729 doublet. 

SW022513 has luminous [OII]A3727, [OIII]AA4959, 
5007, and H/3 line emission, with [OIII] A5007/H/3 = 6. 
[NeV]AA3346,3426 fall outside the atmospheric windows, but 
we detect [Nelll] A3869. We will argue in !|4]that, for galaxies 
with the overall characteristics of SW022513, the [OIII]/H^ 
ratio and the detection of [Nelll] suggest that most of the 
warm ionized gas is photoionized by the AGN. 

Line emission in SW022513 is extended over sizes of 
1.6" X 2.5" in right ascension and declination, respectively. 
The [OIII]A5007 emission-line morphology is shown in the 
upper left panel of Figure [2] The [OIII]A5007 line flux does 
not peak on the continuum peak, but is offset by 0.75" to 
the South (corresponding to a projected distance of about 5 
kpc). 

Emission-line profiles are complex. We identify broad, 
blueshifted components in [OIII] AA4959, 5007 and II/3. 
Unlike [NeIII]A3869 and [Oil] AA3726, 3729, the H/3 and 
[OIII]A5007 lines are fairly bright and do not suffer blending 
of several components with uncertain relative line widths. 
They are therefore particularly suited to investigate their 
line profiles, and we will focus our discussion on those two 
lines. Line profile fits to [OIII]A4959 yield similar results as 
[OIII]A5007 within the measurement uncertainties. 

The broadest [OIII]A5007 component is found near the 
continuum peak, with FWHM=5078 km s"'^ and an offset of 
— 1314 km s~^ relative to the narrow [OIII]A5007 component 
measured in the same spectrum. At the same position we 
also detect a broad component in H/3, with FWHM=1000 
km s~^ and a blueshift of -183 km s~^ relative to the narrow 
component (which has FWHM=369 km s~^). All emission- 
line properties and their uncertainties are listed in Tables [3j 
|4] and [II 



We obtained maps of the relative velocities and line 
widths from fitting emission lines extracted from apertures 
of 3 pixels x3 pixels (0.4" x 0.4") at all positions where 
[OIII]A5007 line emission was detected at >3(t, adopting 
the prodedure outlined in Nesvadba et al. (20081. To ac- 



count for the complex profile of the line, we performed 
fits with 2 Gaussian components for each line. We fit- 
ted the [OIII] AA4959, 5007 doublet simultaneously, requir- 
ing that both lines have the same redshift and line width, 
and a ratio as expected between the two components of 
R(5007,4959)=3. Overall this gives a line fit of 4 Gaussian 
components, where the narrow and broad components of 
[OIII]A4959 and [OIII]A5007 are required to have the same 
kinematics and a line ratio consistent with their Einstein 
coefficients. 

The maps of relative velocities and line widths (full 
width at half maximum) of the narrow [OIII]A5007 com- 
ponent are shown in the upper middle and right-hand panel 
of Figure [2j respectively. Velocities are relatively uniform 
across most of the source, with a sudden redshift jump of 
~200 km s~^ northward of the nucleus. The corresponding 
maps for the broad component are shown in the lower panel 
of Figure [2] We detect a very broad (FWHM~5000 km s"^) 
component of [OIII]AA4959,5007 associated with the con- 
tinuum peak, and a somewhat less extreme, but nonethe- 
less broad component with FWHM~1000 km s~^ associated 
with the peak of the [OIII]AA4959,5007 emission. 

To map line ratios relative to [OIII]A5007 we also fitted 
the spatially resolved emission of H/3 and [Oil] AA3726, 3729, 
where we assumed a single Gaussian distribution for the 
[OH] doublet which is not spectrally resolved. In Figure [5] 
we show the maps of line ratios relative to the narrow 
[OIII]A5007 component. The ratios of [OII]A3727 and of H^ 
to [OIII]A5007 are smallest near the [OIII]A5007 emission- 
line peak, and have larger values near the brightest contin- 



uum emission. 



H/3 and [OII]AA3726,3729 are too faint to map the 
kinematics of multiple components. Therefore we ex- 
tracted integrated spectra associated with the peaks in 
[OIII]AA4959,5007 and continuum emission respectively, and 
from the faint, extended region to the North. These spectra 
are shown in Figure [S] and their properties are summarized 
in Tables |3)|4] and [5] 



3.2.2 SW022550 

SW022550 is a luminous line emitter in the rest-frame UV as 
discussed by [PoUetta et al.| ( |2008| ). Unfortunately at z=3.867 
it is at a somewhat unfavourable redshift for ground-based 
observations, where [OIII]A5007 and [OII]A3726,3729 fall 
outside the atmospheric windows. However, we did iden- 
tify the [OIII]A4959 fine at A=24182 A, corresponding to 
a redshift z — 3.876, which is consistent with the red- 
shifts measured in the rest-frame UV. The [OIII]A4959 line 
IS very broad, FWHM=2212 km s'S not very different 
from the width of the rest-frame UV lines, which have lin e 
widths of up to few 1000 km s"^ ( [Polletta et al.[[2008 |. 
The line ffux is F(4959)=8.2x 10"^^ erg s"' cm ^. For 
F(4959)/F(5007) = 0.3 this corresponds to a [OIII]A5007 
ffux of F(5007)=2.5xl0"^^ erg s"^ cm'^ We did not de- 



tect H/3 but place an upper limit of 4.3 x 10 erg s 
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cm at 3(7 significance, assuming it has a similar width as 
[OIII]A4959. 

We also detect the [NeV] AA3346,3426 doublet, although 
[NeV]A3346 is heavily blended with night sky line residuals, 
making it difficult to measure anything but the line core. 
These lines give redundant kinematic and flux information 
and are well fit with a common redshift of z=3.861 and line 
width of FWHM=2985 km s"^ For [NeV]A3426 we measure 
a flux of 4.3 X 10"^'' erg s"\ and the core of [NeV]A3346 is 
consistent with being 3x fainter as expected from the tran- 
sition probabilities of the two lines. The rest-frame optical 
spectral properties of SW022550 are listed in Table [T] 

All line emission appears compact and associated with 
the compact continuum emission. 



4 TWO QUASARS WITH GIANT 

NARROW-LINE REGIONS AT Z>3.5 

Our spectroscopic results suggest that the AGN is the dom- 
inant source of ionization in our targets. For SW022550 this 
follows directly from the detection of [NeV]A3426. Line emis- 
sion from this galaxy is not spatially resolved, the size of the 
seeing disk of 0.6" implies an upper limit on the radius of 
the NLR of about 2.5 kpc, provided that we did not miss 
fainter, extended structures due to the somewhat unfavor- 



( 2008 1 came to a similar conclusion from the rest-frame UV 



line ratios of SW022550. Roughly solar metallicities have 



able redshift of the source (^3.21. 

At the somewhat lower redshift of SW022513, 
[NeV]A3426 does not fall into the atmospheric windows. We 
do however, observe [NeIII]A3869 and find a relatively high 
[OIII]/H/3 ratio of 6. A priori, both features could be pro- 
duced either by the AGN or by intense star formation in 
galaxies with high ionization parameters and low metalici- 
ties, leading to hot electron temperatures, like, e.g., in high- 



redshift HII or Lyman-break galaxies (e.g., Pettini et al. 


2001 Fosbury et al. 


2003 Nesvadba et al.|2006 2007 20081. 


Villar-Martfn et al. 


2008 1 find that star formation may pho- 



toionize parts of the gas in 3/50 type-2 quasars drawn from 
the SDSS; however, these are much less powerful than the 
targets we observe, and have less luminous optical lines, e.g., 
£{[OIII]) ~ 10"^ erg s"^ compared to fewxlO"* erg s'^ for 
SW022513 (Table|2|. We are not aware of any star-forming 
galaxy without AGN and with £([0/77]) = 10''* erg s"^ 

In addition, the bright 24/im and millimeter continuum 
emission of SW022513 and SW022550 shows these are dusty 
galaxies. Dust and metal lines are very efficient coolants, so 
that for a given ionizing spectrum (e.g., due to intense star 
formation), we may expect lower electron temperatures than 
in HII galaxies and LBGs, which would lead to lower [Nelll] 
fluxes and lower [0III]/H/3 ratios. The FIR/millimeter prop- 
erties of SW022513 and SW022550 are indistinguishable 
from those of submillimeter-selected starburst galaxies at 



z^2, which have supersolar metallicities ( Tecza et al.|2004 
Swinbank et al.|2004| [Nesvadba et al.|2007[) and fall near the 



local mass-metallicity relationship of |Tremonti et al. ( |2004[ ) 
In the absence of an AGN, these galaxies have [OIII]/H/5 ^ 1 
Nesvadba et al.||2007l ), and [Nelll] is 



(Takata et al. 2006 



generally not observed. For SW022550 and SW022513, the 
mass-metallicity relationship would imply a gas-phase oxy- 
gen abundance of 12+log{0 / H) ~ 9.1, supersolar even when 
accounting for the (large) scatter in the relationship (and su- 



also been found by Humphrey et al. ( 2008 1 for dusty, mas- 
sive radio galaxies at high redshift. 

The line emission in SW022513 is well resolved spatially. 
This implies that the QSO narrow-line region (NLR) extends 
to radii of at least 1.25" (corresponding to ~10 kpc). Most 
notably, the peak of the narrow [OIII]AA4959,5007 line is off- 
set by 0.75" (5 kpc) towards South from the continuum peak 
and the peak of broad [OIII] AA4959,5007 emission. A spa- 
tial offset between broad and narrow line emission has pre- 
viously been reported for the submillimeter-selected quasar 
SMMJ1237+6203 at z=2.1 ( [Alexander et al.|2010| |. Without 
additional constraints, this offset could be interpreted as a 
signature of an outflow of turbulent gas, at several 10s of 
kpc from the galaxy traced by the narrow-line emitting gas. 

Our data on SW022513 allow us to perform a more 
complete analysis which suggests a different scenario. First, 
the broad [OIII] line emission is spatially coincident with 
the continuum peak, implying that the emission arises in 
the nuclear region of the galaxy and thus near the AGN. 
We also find that near the peak of narrow [OIII] emission, 
the [Oil] /[OIII] and the H/3/[OIII] hue ratios are smaller 
than near the continuum peak, and hence near the peak of 
broad [OIII] and the likely location of the AGN (Figure 5). 
This implies that the narrow line region is more highly ion- 
ized than the emission immediately surrounding the AGN. 
This increase can be caused by a harder radiation field for a 
constant ionization parameter (ratio of radiation field inten- 
sity to gas density), an increasing ionization parameter, or a 
decrease in gas-phase metallicity as we move from the con- 
tinuum to the narrow-line peak. It could be a combination 
of all three (see Humphrey et al.][2008 for a more detailed 
analysis in related environments). 

Given the bright FIR luminosity of SW022513, the ex- 
tinction in the circum-nuclear region is likely to be high. 
However, large-scale extinction of the emission line gas can- 
not solely determine the variation in [OIII] /[Oil] line ratio. 
The same trend is also found in [OIII]/H/3 and given their 
similar wavelengths (including the [OIII] line at 4959A) the 
influence of extinction must be relatively small. In relation 
to the nebula itself, the effect of dust on the ionization of 
the gas is complex. Since dust provides an additional source 
of cooling and competes directly with the gas for ionizing 
photons due to its large cross section over a wide range of 
wavelengths but also depletes metals from the gas (some of 
which provide important cooling lines like Oxygen and Car- 
bon), a dusty gas will have a different ionization structure 
depending on the gas-to-dust ratio and amount of metal de- 
pletion. So it is not entirely clear what its impact will be, 
but it could well increase the ionization of the gas. 

All these scenarios would produce the situation we ob- 
serve and suggest that the characteristics of the extended 
ionized gas in SW022513 are consistent with a narrow-line 
region surrounding a luminous obscured quasar and extend- 
ing over a radius of ~10 kpc. Radii of few 10s of kpc are 
expected for the narrow-line regions of the most powerful 
quasars ( [Netzer et al.[[2004| ) and consistent with empiri- 



cal constraints for optically selected quasars (Bennert et al 



persolar by 0.4 dex when taken at face value). Polletta et al 



20021 and also type-2 quasars (Greene et al. 20111. In ad- 
dition, the asymmetry of the narrow [OIII] line morphology 
relative to the nucleus supports this scenario as narrow-line 
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regions are often highly elongated and asymmetri c. Given 



the bright FIR luminosity of our targets of 10 Cq ( Polletta 



et al. 2008 1, we expect these are highly dust-enshrouded, 
and that most of the rest-frame UV/optical light is obscured 
on kpc scales, with ionizing and non-ionizing radiation es- 
caping to larger distances only along relatively few lines of 
sight that are comparably free of dust. This is analogous to 
the situation observed for local Seyfert 2s in their extended 
narrow-line regions (e.g., Veilleux et al.||2003 
et al.|2009 for an example at z>2.) 



Lehnert 



This does not imply that the quasar also dominates the 
gas kinematics. Quasar illumination cones without strong 



mechanical effect are known at low rcdshift (Veilleux et al 
2003[ ) as well as at z=2 ( [Lehnert et al.^2009t . Note that this 



applies to radio-quiet quasars. Radio-loud AGN, including 
radio galaxies, often have much more extended emission-line 
regions with very energetic emission-line kinematics (e.g. 
McCarthy et al.||T 996 Villar-Martm ct al. 1999; Baum & 
'McCarthy||2000| 



2006 



2008 



Villar-Martm et al.||2003, ,Nesvadba et al.j 



4.1 Ionized gas mass 

The molecular gas strongly outweights the ionized gas in 
strongly star-forming mergers at z~2 ( Nesvadba et al. 2007 1 , 
but not in powerful radio galaxies at similar redshifts (Nes- 
[vadba et al.|[2008| Nesvadba et al., in prep.) where much 
of the interstellar medium appears photoionized by the 



AGN (e.g., Villar-Martin et al. 1997 Villar-Martm et al 



2003[ [Humphrey et al.||2008[) and accelerated by the radio 



source (e.g., Nesvadba et al.|2008 l. We will now use the H/3 
emission-line flux measured in SW022513 to illustrate that 
the ionized gas mass in this obscured quasar could plausibly 
be as large as the molecular gas mass. 



Similar to Nesvadba et al. ( 2006 2008 1 we will assume 



case B recombination to estimate the ionized gas mass from 
the luminosity of the Balmer lines, setting 



M 



28.2 X 10*" Chpaz n7joo Mq, 



(1) 



where jCh0,43 is the H/3 luminosity in units of 10 erg s 



and is the electron density in units of 100 cm 



This 



relationship is equivalent to Equation (1) of Nesvadba et al, 
( [2006 [ ) assuming a Balmer decrement Ha/H/3 =2.9, and fol- 
lows directly from Oster brock (19891. 



This estimate has two major uncertainties. First, the 
[OII]AA3726,3729 line doublet is blended, so we have no 
direct estimate of the electron density. Second, from H/3 
alone we cannot measure the extinction. We can loosely 
constrain these values with rather generic arguments. First, 
the presence of luminous forbidden emission lines like the 
[OII]AA3726,3729 doublet suggests that electron densities 
are below the critical densities for coUisional deexcitation, 
a few 1000 cm^'^, and we will in the following adopt a fidu- 
cial value of 1000 cm~^ (radio galaxies and submillimeter 
galaxies at similar redshifts have electron densities of few 
100 cm 



Nesvadba et al. 



2006 



2007 



20081 



assume an average extinction of about Ai 



Second, we will 
^ 2 mag, which 



Hq:/H/3 line ratios in the extended ionized gas of z~2 ra- 
dio galaxies. For SW022513, this gives an H/3 luminosity of 
Chi3 = 7 ~ xlO*'* erg s~^, a factor ~ 10 higher than the 
observed value (Table [2]|. 

With these two assumptions, we find an ionized gas 
mass of 2 X 10^° M© (compared to 2 x 10^ M© if we strictly 
use the observed H/3 flux). This can obviously only be an 
order-of mass estimate, but illustrates that the mass of ion- 
ized gas in SW022513 is likely between 10% and 100% of 
the molecular gas mass. This is in the same range as found 
for radio galaxies, and a factor 100—1000 higher than for 
strongly star-forming galaxies. For SW022550 we only have 
an upper limit on the H/3 flux, which is not very constrain- 
ing. The above considerations suggest that the ionized gas 
mass estimate and derived quantities can only be accurate 
at an order-of-magnitude level. 



4.2 What is the systemic redshift? 

The lack of robust measurements of the systemic redshift 
is a major uncertainty in all studies of the gas dynamics 
in galaxies at high redshift, and yet, is indispensible to 
interpret the emission-line kinematics. At distances where 
direct measurements of stellar absorption line kinematics 
in AGN hosts are possible, observations suggest that the 
narrow emission-line components show only moderate off- 
sets from the systemic velocity. This includes Seyfert galax- 



ies (Nelson & Whittle 19961 as well as powerful, type-2 
quasars in the SDSS ( Greene fc Ho|2005a \ , and also galaxies 
with very broad blueshifted [OIII]AA4959,5007 components 
( Wong et al.|2006[ ). 

The narrow [OIII]A5007 component in SW022513 has a 
relatively small velocity gradient, Av ~200 km s^^ (Fig- 
ure |2|, and lines of different species of ionized gas have 
very similar velocities. This is broadly consistent with what 
would be expected for rotation or merger-driven kinemat- 
ics of galaxies with stellar masses of ~ 1O"M0 and IS m 



the velocity range of submillimeter galaxies at z^2 (S win- 



bank et al. 2006 Nesvadba et al. 20071 as well as nearby 
ULIRGs ( [Cohna et al.[[2005[ ). We wiU thus in the following 
assume that the narrow optical emission lines of SW022513 
trace the systemic redshift within ^ 100 km s~^. Specifically 
we will use the redshift of the narrow H/3 component near 
the continuum peak (which we identify as the nucleus, sQ, 
z=3.4247 (Table [sjl to approximate the systemic redshift. 
H/? is measured at good signal-to-noise (SNR=8.5(t) and 
is more representative for the overall gas kinematics than 
[OIII]AA4959,5007, which are very sensitive to ionization ef- 
fects. 



4.3 Blue wings as signatures of outflows 

If the narrow lines of AGN host galaxies are approximate 
tracers of the systemic velocity, then the broad, blueshifted 
components often found in [OIII]AA4959,5007 as well as 
other lines, most likely trace gas that is in outflow (e.g.. 



is the average in submillimeter-selected galaxies without 


Heckman et al.||1981| [Greene & Ho[[2005b 


[Morganti et al. 


strong AGN (Small et al.|2004l that are the closest analogs 


2005 


Nesvadba et al.|2007 Komossa et al. 


2008 


Holt et al. 


to our galaxies in the far-infrared/submillimeter (and con- 


2008 


INesvadba et al. 2008| [Greene et al 


1 2009 


[Spoon & 


sistent with a nuclear Av ~4.6 mag derived by|Polletta et al.| 


Holt 


2009; Spoon et al.||2009[ [Nesvadba et al.||2010l. Scal- 



2008). Nesvadba et al. (20081 find Av =1-4 mag from the ings between the [OIII1A5007 width and AGN power (either 
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radio power; |Hecknian et al.|198l| or bolometric luminosity; 



Greene & Ho 



2005b I broadly support this picture. Follow- 



ing these previous analyses, we therefore consider the broad 
wings of [OIII] and other lines in SW022513 and the broad 
[OIII] lines in SW022550 as signatures of outflows. 

In SW022513 the broadest hue widths of FWHM^5000 
km s^^ seen in [OIII] are broader than the other lines in the 
same aperture, and broader than all other lines including 
[OIII] in all other apertures. The [OIII] emissivity is very 
sensitive to electron temperature and excitation conditions, 
but not the total gas mass, so this component is unlikely to 
trace large amounts of gas. The most strongly blueshifted 
gas seen in the more representative H/3 line, blueshifted by 
of order —1000 km s"'^ relative to systemic, likely gives 
a more robust estimate of the overall gas kinematics (see 
also i4.4l. We ran a Monte-Carlo simulation to infer for 



rate to about 20 km s~^. In the millimeter, wavelengths are 
measured relative to a local oscillator, therefore uncertain- 
ties in wavelength calibration are much smaller and can be 
neglected. 

Formally, our Gaussian fits imply a velocity offset of 
-183±67 km s^^ for an assumed Gaussian line core of the 
blueshifted H/3 component, and of -181±47 km s'^ for the 
Gaussian core of the CO(4-3) line. Molecular line emission 
is therefore found with blueshifts of up to —900 to —1000 
km s~^ relative to the narrow H/3 component, which, as 
argued in §4.2| is our most robust measure of the systemic 
redshift. Regardless of the detailed match between molecular 
and ionized gas, and neglecting possible projection effects, 
this is likely to be more than the escape velocity (j|5.3[). 



what H/3/ [OIII] A5007 line ratio we could have detected a 
broad H/3 component with FWHM=5000 km s"^ and at 
the same redshift as the broad [OIII]AA4959,5007 lines (see 
§ 3.1 of [Nesvadba et al.||2011| for more details). For 1000 
throws we found that H/3 would have been detected at 3a for 
[OIII]A5007/ up s; 6.6, compared to [OHI] A5007/H^=3.8 
for the narrow component in the same spectrum. If the 
assumptions of §4.1| hold, then this would imply a 3a up- 
per limit to the mass of high-velocity ionized gas (with 
FWHM=5000 km s"^) of 1.5 x 10** Mq (neglecting extinc- 



5 THE WIND IN A Z=3.5 QUASAR 

Our targets have all the hallmarks of being in a short, deci- 
sive, and very complex stage of their evolution, where star 
formation, AGN, and gravitational interactions are likely 
to release an energy of IQ^^"^" erg, the equivalent of the 
binding energy of the host galaxy and its dark-matter halo. 
At their current luminosities, the AGN in SW022513 and 
SW022550 may release such energies in a few lO''"^ yrs. 



tion), or 1.5 x 10^ Mq (for Ay =2 mag; see i|4T|, roughly which corresponds to the lifetimes over which A GN may 



10% of the mass found at more moderate velocities. Inter 
estingly, the higher [OIII] /H/3 ratio suggests this gas is more 
highly ionized than the narrow- line emitting gas. 

Notwithstanding, H/3 also shows the most strongly 
blueshifted components near the nucleus (see Tables |3] to 
[5|. This is expected if the gas flows are driven by the AGN 
and decelerate as they interact with ambient gas at larger 
radii. 



4.4 Comparison with molecular gas 



maintain bolometric luminosities > lO^'^ erg s ^ (|Hopkins 



et al. ] |2005| ). Similarly, the canonical energy release of 10"^ 



PoUetta et al. (20111 recently discussed integrated C0(4- 



3) emission-line spectroscopy of SW022513 and SW022550 
obtained with the IRAM Plateau de Bure Inteferometer at 
a spatial resolution of ^4". Both galaxies have luminous 
millimeter CO(4-3) line emission that corresponds to --^ 4 x 
10^" Mq of molecular gas, assuming a factor of 0.8 Mq K km 
s~^ pc^ to convert CO luminosity to a molecular gas mass. 
SW022550 has a double peaked profile with two components 
separated by 440 km in velocity. SW022513 has a single, 
very broad component with FWHM~1000 km s~^, which is 
not very well fitted with a single Gaussian. 

In Figure |6] we compare the CO (4-3) and [OIII]A4959 
line profiles of SW022550. The redshift of [OIII], 
z=3.876±0.001 falls between that of the two CO peaks. 
In SW022513, the CO(4-3) line profile shows an excellent 
match with the H/3 profile extracted from the nuclear aper- 
ture centered on the continuum peak (Figure [7|. This in- 
cludes in particular the broad, blueshifted wings with veloc- 
ities of up to -1000 km s~^ relative to the narrow H/3 compo- 
nent. As we argued in j ]4.2| the narrow H/3 component most 
likely approximates the systemic velocity to about 100 km 
s-i. 

Comparison with telluric night-sky lines suggests that 
the absolute wavelength scale in our SINFONI data is accu- 



erg per supernova explosion corresponds to a total energy 
of about xlO^* erg released in the formation of 10^^ Mq 
of stars, and observations have shown that the hydrostatic 
pressure of intense starbursts may balance gravity in z=2 
galaxies very similar to our targets ( [Nesvadba et al.||2007[ ). 
Major mergers naturally release the equivalent of the bind- 
ing energy of the merger remnant during the interaction. 
Hence, each of these mechanisms could have a strong influ- 
ence on the kinematics and thermal state of the gas in our 
targets. What mechanism is dominating the gas dynamics 
in SW022550 and SW022513? 

Important to address this question is not only the to- 
tal energy injection, but also the timescales over which this 
energy is released, and the efficiency with which the energy 
output is turned into an input of (thermal or mechanical) 
energy into the ambient gas. Narayanan et al. (2008j) find in 
SPH simulations of a gas-rich, major merger associated with 
intense star formation and AGN activity, that the AGN af- 
fects the gas more drastically than star formation and grav- 
itational collapse because of the shorter energy injection 
time. Star formation and gravitational interaction release 
their energy in few 10* yrs and up to 10^ yrs, respectively, 
much longer than the 10®~^ yrs lifetime of the AGN. 

Observations lead to a similar conclusion. The spectral 
signatures of winds in submillimeter galaxies without strong 
AGN component appear much more subtle, with smaller 



blueshifts, and a smaller mass of entrained material (Nes- 
'vad ba et al.|2007 l. Overall, these galaxies have velocity gra- 
dients and line widths of up to a few 100 km s~^, and more 
regular line profiles (e.g., Swinbank et al.|2005 2006 Takata 



et al. 2006 1 in spite of stellar masses and star-formation 



rates very similar to our targets. Low-redshift ULIRGs in- 
cluding very advanced mergers have velocity gradients and 
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line widths in a similar ranee ( [Colina et aTlpOOSt , and are 
much smaller than ~1000 km s~ at any stage of merging. 

SW022513 and SW022550 both show signatures that 
the AGN does afltect the warm ionized ISM, evidenced 
through the broad, blue components in SW022513 and 
the very broad FWHM of [OIII]A4959 in SW022550 
(FWHM>2000 km s'^). However, FWHM>1000 km s'^ are 
only found around the nucleus in both cases, and , since we 
cannot constrain extinction very well, it is not clear if the 
warm ionized gas represents a major fraction of the ISM 



lono et al. 2007 Feruglio et al. 2010 Irwin et al. 20101. These 



(f 4.1 1. The CO profiles of both galaxies, which trace about 
4 X 10^" Mq in molecular gas in each source, are very dif- 
ferent in SW022513 and SW022550 ( [Polletta et "ar]|20TT| ). 
In SW022550 the CO profile resembles the fairly common 
double-peaked profiles found, e.g, in submillimeter galaxies, 
where they are commonly interpreted as signatures of merg- 
ers or rotating disks (and where even disks may be signatures 
of (advanced) mergers; pownes fc Solomon||1998|) m agree- 
ment with the merger models of Narayanan et al. ( 2008 1. In 



either case, the barycenter falls roughly inbetween the two 
CO peaks. This could suggest that the AGN affects only a 
small part of the multiphase ISM in SW022550, traced by 
the broad [OIII] hue. 

SW022513 is however different. The CO line profile is 
irregular and neither a clear single nor double-peaked Gaus- 
sian. The line extends to large rela tive velocities of up to 
— 1000 km s~^ from systemic {{4.4\, similar to the ionized 
gas near the AGN (Fig. [7|. Outflows driven by the AGN 
are the most common interpretation of blue, broad wings 



of ionized gas (S4.3I, which may suggest that parts of the 
molecular gas in SW022513 may be tracing outflowing gas 
as well. Winds are inherently multi-phase phenomena, where 
the warm ionized gas seen through optical line emission is 
being entrained by a hot, tenuous medium. It may there- 
fore a priori not be entirely surprising to find an associated 
phase of molecular gas. However, this does raise fundamental 
questions of how molecular gas, which is 1-2 orders of magni- 
tude denser (iV ~ 10^~* cm"^) than ionized gas (TV ~ 10^"^ 
cm~^) is accelerated to high velocities. This could be allevi- 
ated if much of the gas is diffuse and distributed across the 
galaxy, perhaps through tidal effects or a starburst-driven 
wind ( Narayanan et al.|2008 l. 

It is also possible that the molecular gas is forming m 
situ in the outflow, as suggested by recent studies of the de- 
tailed mass and energy exchange in turbulent, multiphase 
gas in nearby extragalactic environments with galaxy-wide 
shocks, including radio galaxies ( [Guillard et al.|[2009 Nes- 



Ivadba et al. 2010, Ogle ct al. 2010] see also Krause & Alexan- 
der, 2007 ). In these scenarios, the ionized line emission 
may arise from the turbulent mixing interfaces associated 
with the same clouds, which would explain why ionized and 
molecular gas have similar velocities. |Papadopoulos et al.| 
| 2010l [2008| ) recently found bright CO(6-5) and CO(3-2) 
line emission in the nearby radio galaxy 3C293, which is 
a posterchild of a jet-driven outflow ( Morganti et al.||2005 



Emonts et al.|p"005 l. CO line ratios in 3C293 suggest that 



most of the molecular gas is dense, turbulent, and gravita- 
tionally unbound, as expected in a multiphase scenario. 

Possible kinematic signatures of outflows from star- 
forming AGN host galaxies traced through CO line emission 
have previously been reported in a few nearby cases (|Apple-| 
|ton et al.|[2002l [Sakamoto eFaL][2006l [Alatalo et al.||2010[ 



galaxies have AGN power and star-formation rates that are 
lower by 1-2 orders of magnitude than in SW022513. Outflow 
velocities are typically a few 100 km s~^, and entrainment 
rates can be up to few 100 Mq yr~^. The CO profiles found 



by Sakamoto et al. (20061; Alatalo et al. (20101 are weU 



matched by HI absorption-line profiles tracing neutral out- 
fiows, and providing robust evidence for the outfiow interpre- 
tation. However, typically these winds include only a small 
fraction of the CO luminosity, unlike in SW022513 where 
most of the CO emission appears to be blueshifted. This 
may partially be due to the higher-excitation gas probed in 
the J=4-3 transition, or it may imply that more gas is being 
entrained compared to low-redshift galaxies. This is not im- 
plausible, given the much greater AGN power in SW022513 
compared to nearby Seyfert galaxies, and the large molecu- 
lar gas mass. We will in the following quantify the necessary 
and available amounts of kinetic energy in SW022513 to in- 
vestigate whether the AGN may plausibly drive the gas to 
the velocities observed. 



5.1 Kinetic energy of the gas 

In Figure [T] we compared the line profiles of H/3 and CO, 
and argued that, by analogy with a large number of previous 
AGN studies, the systemic velocity is typically well approx- 
imated by the narrow component of optical emission lines, 
whereas the outfiowing component is in the blue wing. To 
estimate the kinetic energy in the outflowing gas, we there- 
fore decompose the CO and H/3 line profiles into a 'systemic' 
and a blueshifted component, finding that about 40% of the 
total CO (4-3) emission-line luminosity is in the outflowing 
component (and about 30% of H/3). Assuming that the CO- 
to-H2 conversion factor and the CO excitation are indepen- 
dent of velocity, this suggests that 40% of the molecular gas 
is in the blueshifted component. 

Based on these assumptions, we estimate the kinetic 
energy directly from the line proflle, by setting 



Eh 



1/2 M, 



(2) 



where Mi is the gas mass in each velocity bin Vi (relative 



to the systemic velocity, S4.2I, and including only the flux 



from the blueshifted component. This corresponds to a ki- 
netic energy of about 4 x 10^* erg. Adding the ionized gas 
mass would add another 0.2 — 2 x 10^* erg, depending on 



extinction (see S4.ll. Obviously, these estimates have large 
systematic uncertainties related to the molecular gas mass 
estimate, excitation conditions, and velocity estimates in- 
cluding projection effects. We therefore consider this esti- 
mate accurate to about an order of magnitude. 



5.2 Is the AGN capable of driving the outflow? 

At the beginning of the section we have argued that a 
starburst-driven wind and an interaction are unlikely to ac- 
celerate signiflcant amounts of gas to velocities of —1000 km 
s~^. We will now investigate for two popular AGN feedback 
mechanisms, if they may plausibly explain the observed gas 
kinematics. 
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5.2.1 Radiation pressure from AGN and starburst 

Radiation pressure is often invoked to explain how AGN 
may expel significant amounts of gas from their host galaxies 
(e.g., |King|[2003l [Murray et aL]|2005[ ). [Murray et al.| ( |2005 l 
discussed the outflow velocities that can be produced by 
radiation pressure from AGN and intense starbursts, finding 
that (their Equation 17) 



V{r) = 2aJ{ 



C 

Cm 



l)ln 



Ro 



(3) 



where a is the stellar velocity dispersion of the host galaxy, 
and C is the quasar luminosity. Rq is the launch radius of 
the outflow, and r the radius at which the velocity of the 



wind is measured. Ca 



4 /s 



is a critical luminosity 



that depends on the stellar velocity dispersion cr, the speed 
of light, c, gravitational constant, g, and the gas fraction, 
fg. For C > Cm, radiation pressure may launch a wind. 
These equations are appropriate for the limiting case of an 
optically thick wind, in which case the interaction is most 
efflcient. 



PoUetta et al. ( 2008 1 estimated the stellar mass and 
luminosity of the starburst and AGN from the dust and 
stellar emission of SW022513 with an exquisite set of 
multi-waveband photometry. They find a stellar mass of 
MsteHar ~ 2 — 4 X 10^^ M©, aud bolometric luminosities 
of 6 X 10*® er(;s~^ and 4.8 x lO'"^ erg s~^ for the starburst 
and AGN, respectively. For pressure-supported galaxies with 
approximately isothermal mass profile, this mass range cor- 
responds to stellar velocity dispersions of cr = 300 — 350 km 
s~^. This can be found from setting M= ccr^ r^/G, where 
a is the stellar velocity dispersion, re the effective radius, 
M the stellar mass, and G the gravitational constant. For 
our calculations we adopted re ~ 2-3 kpc, and c = 5 (Ben 



der et al.nl992 ). To give a lower limit on the gas fraction, we 



use the molecular gas mass estimate of [Polletta et al.[ ( [2011[ ) , 
4x10^° M0, which gives a gas fraction of order fg = 0.1 — 0.2 
for a stellar mass of Msteiia 
[et al.|2008| |. 



= 2 - 4 X 10'' Mq (PoUetta 



Following Murray et al. (20051, to launch a wind in 

4 X 10" 



a galaxy with Msteiiar 



350 km s 



and M„ 



2 - 4 X 10'' Mq with cr=300- 
= 4 X 10^° Mq would require 
Cm = 2.5 — 3 X lO*'^ erg s~^. This is the most optimistic case 
consistent with our observational constraints, and would re- 
quire a bolometric luminosity that is about 2.5— 3 x greater 
than observed. To accelerate a wind to a terminal velocity 
of 1000 km s~^ (as suggested by the broad blueshifted 
components in SW022513) would then require a bolometric 
luminosity of at least 4.5 — 5 x 10*^ erg s~^ for a galaxy with 
a =300—350 km s~^, a factor 4-5 larger than what is mea- 
sured. Thus, the luminosity of SW022513 is at least a factor 
4 lower than required, if we consistently use the most opti- 
mistic estimates implied by our observations. It can be up 
to factors 10 lower for other plausible choices of parameters. 

For these estimates we assume a launch radius of the 
wind, -Ro, of a few 100 pc (the sizes of the circumnuclear 
molecular disks found in low-redshift ULIRGs IDownes fcl 
|Solomon||1998[ and th e lowest value in the AGN feedback 
models of Ciotti et al.[[2009 20101, and an outfiow radius, 
r, of 5 kpc, which roughly corresponds to the radius of 
the narrow-line region in SW022513. A larger launch ra- 
dius (perhaps more plausible given the large gas masses) or 



a larger size of the outflow region, enhances the required en- 
ergy. This suggests that radiation pressure may only drive 
the outflow if all parameters take values that are strictly 
at their lower limits and if we underestimate the luminosity 
by at least a factor 4. Otherwise, driving an outflow like in 
SW022513 through radiation pressure as proposed by |Mur-] 
ray et al. ( 2005 1 would require an AGN and starburst with 



at least factors of a few higher bolometric luminosity than 
in SW022513, although this is one of the most luminous ob- 



scured high-z AGN in the SWIRE survey (PoUetta et al 



20081 



5.2.2 Radio source 

Examples of outflows in the literature associated with radio- 
loud AGN are numerous, and include not only powerful radio 
sources (e.g.,|Morganti et al. 12005 Holt ot al. 2008, ^Nesvadba 



et al.||2008| |Spoon & Holt 


2009; ISpoon et al.[[2009| [Nes- 


vadba et al. |2010l, but also Seyfert galaxies with low radio 


power (Capetti et al. 1999 


Ulvestad et al. 1999 Reynolds 


et al. 2009 Barbosa et al. 


2009). Weak radio sources are 



always easily identified, and typically spatially asscociated 
with non-circular gas motions. The same is found for radio- 



quiet quasars (e.g., Leipski & Bennert 2006 Leipski et al 



2006 1 . Note that "radio-quiet" does not mean "radio-silent" , 
but that the radio power is less than 10 x greater than the 
optical luminosity of a quasar, irrespective of whether the 
optical continuum is dominated by the direct or scattered 
AGN light, nebular continuum emission, or the obscured or 
unobscured stellar continuum of the host galaxy. 

Even relatively low-power radio sources can acceler- 
ate gas to high velocities of ^1000 km s~^ (e.g. Capetti 
et al.[[1999[) a nd produce equally large line widths (e.g., 
Capetti et"ar[[1999[ [Reynolds et al.[[2CI09[ ). [Heckman et al. 
U981[ ) fouiidthat radio )ower correlates with the FWHM 
of the [OIII]A5007 line. Measuring the [OIII]A5007 FWHM 
in SW022513 from the integrated spectrum, we find that 
it falls well within the scatter of the IHeckman et al.l rela- 
tionship (Figure [sl . The same is true for SW022550 and 



SMM J1237+6203 (with the FWHM of [Alexander etal 
20T0| 



Weak radio sources where the internal pressure does not 
greatly exceed the ambient pressure, may deposit most of 
their energy in the ambient medium ( Gallimore et al.|2006 \ 
and may thus be more efficient in infiating bubbles than 
more powerful radio sources, where only a few percent of the 
kinetic jet energy is being deposited in the ambient medium 
(e.g., [Nesvadba et al.|2008l ). [PoUetta et al.| ( [2008[ ) estimated 
the radio power of SW022513 (and SW022550) to be of order 
10^^ W Hz"^ This roughly corresponds to the FRI/FRII 
divide at low redshift, and appears as moderately strong only 
if compared to the very powerful radio-loud quasars and 
radio galaxies at high redshift. At more moderate redshifts, 
radio sources with fairly similar power can trigger significant 
outflows (e.g. [Morganti et al.|2005} [Spoon fc Holt|2009[ ) and 
have a profound impact on the molecular gas in their hosts 



(e.g.,jPapadopoulos et "ar|[2010[ [Nesvadba et aL]|2010[ [Ogie 



et al.|2010|[Alatalo et al.|2010[ ). 

The interaction efficiency of a radio source with the am- 
bient gas depends critically on the gas density (e.g., |Capetti[ 
et al.|[l999 l, which in high-redshift galaxies is likely higher 
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than at low redshift. For example, De Young ( 1993 1 suggest 
that trapping a radio source with about 10'^'^ W Hz~^ in 
the ambient gas requires a volume-averaged gas density of 
order 100 cm~^. Assuming a volume-averaged density may 
appear somewhat artificial, but has previously shown to pro- 
vide reasonable constraints ( De Young||1993 Capetti et al 



1999 1, and is certaintly well matched to the crudeness of our 



observational constraints. For low-redshift galaxies, plausi- 
ble estimates of the vol ume-averaged am bient gas density 
are of order 1-10 cm^'^ (De Young 19931, however, a value 
of 100 cm"'^ is very similar to the average density 50-380 
cm"^ in SW025513 and SW022550 implied by the molecu- 
lar mass of 4 X 10^'' Mq ( [Pollettaet al.|201l[ ), and assummg 
a spherical gas distribution with radius 1-2 kpc and a filling 
factor of unity. As a consequence, the higher gas densities in 
high redshift galaxies could boost the effect of moderately 
strong radio sources on their surrounding gas compared to 
nearby galaxies, where such jets escape more easily. 

In many similarly strong, low-redshift radio sources the 
total gas mass involved in the outflow is much smaller than 
in SW022513 (e.g. 



Morganti et al. 



2005). The H/3 lumi- 
1000 km s"^ im- 
^ Mo, estimated 



nosity in the component with FWHM= 
plies an ionized gas mass of 1.2 X 10** ivi0 
with the method, assumptions, and uncertainties presented 
in §4.1[ We therefore test explicitly if the radio source 
in SW022513 may provide enough mechanical energy to 
produce the observed emission-line kinematics of H/3 and 



CO. Willott et al. ( 1999 1 estimate the jet kinetic energy 
from the observed radio power £151,28 at 151 MHz, setting 
Cme^h = 3 X 10^* /w^-^isMS W. £151, 28 is the observed ra- 
dio power at 151 MHz in units of 10^® W Hz"^ sr~^, and fw 
is a fudge factor taking into account the most salient astro- 
physical uncertainties. Typically fw 



10 (see also |Cattaneo| 
et al.|2009 l. Estimating the 151 MHz radio luminosity from 



the measured radio fluxes at 1.4 GHz and 610 MHz (PoUetta 
|et al.|[2008| , we find a mechanical energy injection rate of 
2 X lO^^^erg s'V 

To estimate the total kinetic energy released by the ra- 
dio source, we need to estimate the AGN lifetime. |Blun-| 
dell fc RawIIngsl ([T999^ argue that the most powerful ra- 
dio sources at high redshift must be very young, about 10^ 
yrs. We will use this estimate as a lower limit on the age 



of the radio source. However, Sajina et al. (20071 find that 



1/3—1/4 of intense infrared-selected starbursts at zjs2 have 
moderately radio- loud AGN, with radio power and star- 
formation properties broadly similar to our sources. This 
suggests that moderately powerful radio sources are not un- 
common amongst intensely star forming galaxies at high red- 
shift. Assuming that all such galaxies are moderately bright 
radio sources at some time of their evolution yields a lower 
limit on the lifetime of radio sources in individual galaxies 
relative to the lifetime of the starburst. For typical star- 
formation timescales of few 10* yrs, this would suggest that 
moderately radio-loud AGN at z~2 have longer radio life- 
times than the most extreme radio sources, up to of-order 
10* yrs, more resembling the long activity periods of weak 



radio sources in nearby Seyfert galaxies (Gallimore et al. 
2006 1 than the very short lifetimes (of-order 10' yrs) of very 



powerful radio galaxies at high-z (e.g., Blundell & Rawlings 



1999) 



erg s~^ in SW022513 were only maintained for 1 x 10^ yrs, 
the radio source would provide sufficient mechanical energy, 
7 X 10^* erg, to accelerate large amounts of of molecular 
and ionized gas as is observed, if all of the jet's mechanical 
energy is transferred to the gas at a conversion efficiency 
of 100%. For longer timescales, the same could be achieved 
with lower conversion efficiencies, of-order few 10%, closer 



to those found in very powerful radio galaxies at z^2 (Nes- 
vadba et al.|2008| ). 



5.3 Will the gas escape? 

AGN-driven outfiows were postulated by galaxy evolution 
models to sweep up and unbind the remaining gas at the end 
of merger-triggered major episodes of star formation at high 
redshift. Major gas-rich mergers are the most plausible pro- 
cesses that may trigger intense star formation and AGN ac- 
tivity in high-redshift galaxies, however, our data do not re- 



veal multiple continuum components for either source (S 3.1 ) 



This could imply that both components are at distances less 
than the spatial resolution of our data (~4— 5 kpc) or that 
their stellar components have low surface brightnesses, per- 
haps because they are being tidally disrupted. Either in- 
terpretation would be consistent with an advanced merger 
stage. The same is suggested by the short gas consumption 
timescales of few x 10^ yrs, about 10% of the age of the 



starburst (Pollettaet al. 20111 



Finding outflows of molecular and ionized gas associ- 
ated with an obscured AGN is therefore certainly consis- 
tent with this broad evolutionary picture. However, a crit- 
ical question is whether the gas will escape. In §5.2. 1| we 
estimated that the velocity dispersion corresponding to the 



stellar mass of SW022513 derived by |Polletta et al.|p008 l 
is likely about a =300-350 km s~^, which would imply an 
escape velocity \/2 a ~ 500 km s^^. 

Finding gas with velocities of up to 1000 km s"'^ cer- 
tainly suggests that at least a fraction of the outfiow in gas 
may ultimately become unbound, depending on how much 
energy is being dissipated by the outflow before the gas has 
reached large radii. Recent observations of molecular gas 
in nearby radio galaxies finding roughly similar amounts of 
thermal and kinetic energy suggest that dissipational pro- 
cesses in the molecular gas could be important for the dy- 
namics of molecular winds ( Nesvadba et al.|2010 \ . The most 
direct evidence for gas escaping from the galaxy would there- 
fore be the detection of outfiowing gas that extends well 
beyond the radius of the host galaxies, with large velocity 
offsets, and aligned with the radio axis, as has been found 



in powerful radio galaxies at z^2 (Nesvadba et al. 2006 



Even if the current energy injection rate of 2 x 10 



|2007[|2008]|2011[ ). At any rate, our observations suggest that 
"radio-quiet" , powerful AGN at high redshift may accelerate 
significant fractions of their interstellar medium to velocities 
near or above the escape velocity even if their radio sources 
are rather unconspicuous. 



6 SUMMARY 

We present an analysis of deep rest-frame optical integral- 
field spectroscopy of two powerful obscured quasars at 
z^3.5, SW022513 and SW022550. These are the most lumi- 
nous 24^m emitters in the SWIRE survey at Js2 and have 
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previously been discussed by Polletta et al. (2008 20111 



Our main results are as follows: 

(1) The optical line emission in both galaxies is dom- 
inated by the luminous narrow-line region ionized by the 
hard quasar spectrum. Emission lines are very broad in 
SW022550, FWHM=2200 km s'^ for [OIII]A4959. Emission- 
line profiles in SW022513 are complex. For example, 
[OIII]A5007 is dominated by a 'narrow' component with 
FWHM=1275 km s~^ and has a broad, blueshifted com- 
ponent of up to FWHM=5000 km s"^ These hue widths 
reflect the kinematics of the narrow-line region, the broad- 
line region is obscured. 

(2) In SW022513 the line emission is spatially extended. 
We identify the nucleus with the continuum peak and site of 
broadest [OIII]AA4959,5007 line emission (FWHM=50Q0 km 
s~^). The peak in narrow [OIII]AA4959,5007 line emission is 
offset by 0.75" to the South (~5 kpc). This suggests that the 
ionized gas reaches the largest velocities near the nucleus, 
and is surrounded by an extended narrow-line region. For 
bright AGN like in SW022513 narrow-line region sizes of 
10s of kpc are possible. 

(3) For S W022513, and compar ing with the CO (4-3) 
observations of |Polletta et al.] ( |2011[ ), we find that the ion- 
ized gas mass amounts to 10%-100% of the molecular mass. 
CO (4-3) and H/3 line profiles are well matched suggesting 
both may originate from the same gas clouds. Using the 
narrow H/3 component near the nucleus to define the sys- 
temic redshift (as commonly done in AGN host galaxies) 
we find that about 40% of the CO line emission is from 
gas that does not participate in systemic motion, but is 
blueshifted to velocities of up to —1000 km s~^. Such large 
velocity offsets are not suggested by merger models including 
starburst-driven winds and direct empirical evidence, but 
could be a signature of AGN-driven winds as postulated, 
e.g. by [Narayanan et al.|p008l l. In SW022550 the CO(4-3) 
line profile is double-peaked and very different from that of 
[OIII], which may suggest that the bulk of molecular gas 
is not affected by the AGN in a similar way, although the 
FWHM of the [OIII] hue of >2000 km s~^ suggests that the 
AGN interacts with the ISM nonetheless. 

(4) SW022513 and SW022550 host AGN with bolomet- 
ric luminosities of ~ 5 x 10*® erg s~^ and moderately strong 
radio sources. Comparing with the expected characteristics 



of radiation-pressure driven winds (following Murray et al 



2005 1 and mechanical energy injection from the radio source 
(following [Willott et"ar]|1999|) we find that it is difficult 



to produce the observed velocities with radiation pressure, 
whereas observations and basic energy considerations sug- 
gest the radio source as a possible driver. In this case, mod- 
erately radio-loud ULIRGs like SW022513 and SW022550 
could be 'scaled-up' versions of nearby Seyfert galaxies and 
'radio-quiet' (but not radio-silent) quasars with weak radio 
sources. 
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Figure 1. Integrated spectrum of SW022513 showing H/J and the [OIII] AA4959,5007 doublet. Note the broad [OIII] components. Hatched 
hnes show Gaussian fits with parameters given in Table [2] and are shifted by an arbitrary value along the ordinate. The lower panel 
shows the fit residuals. The inset shows [OIIJAA3726, 3729 (the doublet is blended due to the broad line widths and is fitted with a single 
Gaussian profile for simplicity) and [NeIII]A3869 in the H-band part of the spectrum. 
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Figure 2. [OIII]A5007 emission-line maps for SW022513 for the narrow (top panel) and broad component (bottom panel). From left to 
right, we show the maps of integrated fluxes, velocities relative to the velocity at the positio of the continuum peak, and line widths. 
Color bars indicate the velocities and line widths in km s~^. Contours show the continuum morphology in all panels. 
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Figure 3. Spectra of SW022513 extracted from 0.4" X 0.4" box apertures centered on the continuum pealc (upper panel), the peak in 
[OIII]A5007 emission- line surface brightness (mid panel). The lower panel shows the integrated spectrum of the faint emission-line region 
north of the nucleus (extracted from an aperture of 0.4" x 0.5" in right ascension and declination, respectively). 
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Figure 4. Integrated spectrum of SW022550 showing [NeV]AA3345,3426 (top panel) and [OIII]A4959 (bottom panel). [OIII]A5007 falls 
outside of the K-band atmospheric window, and H/3 is not detected. Hatched lines mark Gaussian line fits, and the bottom panels show 
the fit residuals. [NcV]A3345 is superimposed by several night sky-line residuals, and we did not fit the line, but plot the line fit with the 
redshift, and line width expected from those measured for the [NeV]A3426 line, and use a line ratio of R(3426,3345)=2.8 expected from 
the transition probabilities of each component of the doublet. 
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Figure 5. Line ratios in SW022513 (left to riglit): Flux ratio of the broad to narrow component of [OIII]A5007. Ratio of HI3 to [OIII] A5007 
(narrow component). Ratio of [Oil] AA3726, 3729 to [OIII]A5007 (narrow component). Both for H/3 and [Oil] line ratios relative to [OIII] 
decrease with increasing distance from the AGN. Contours show the continuum morphology in all panels, corresponding to the position 
of the AGN. 
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Figure 6. Comparison of the line profile of CO(4-3) (black line'PoUetta et al.'2011') and [OIII] A4959 (blue line) in SW022550. The lower 
panel shows the residual after fitting the CO profile with two Gaussian distributions (which are shown in red in the upper panel). The 
[OIII]A4959 profile was scaled arbitrarily. 
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Figure 7. Comparison of the line profile of H/3 near the nucleus (upper panel, see also the mid panel of Figure [sjl and the CO(4-3) 
emission-line spectrum of SW022513 (lower panel) presented by |Polletta et al^ l |2011[ |. The red, cyan, and blue Gaussian distributions show 
the CO(4-3), narrow H/3 component, and sum of both spectra, respectively. The red and cyan distributions have the same wavelength 
and width in both panels. 
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Figure 8. Radio power measured at 1.4 GHz in the rest-frame as a function of FWHM of tlie [OIII]A5007 emission line. Small black 
dots show nearby AGN spanning a large range in radio power, and are taken from the original version of this plot shown in|Heckman| 
|et al.| | |198l| |. Large blue dots show SW022550, and SW022513, the red triangle shows SMMJ1237+6203 ( [Alexander et al.pOlof 
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Line ID 


Ao 


z 




FWHM 


FWHM 


flux 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


[OIII] 


4959 


3.876±0.001 


24182±7 


179.5 


2212±178 


8.2±0.76 


[NeV] 


3456 


3.861±0.002 


16653±6 


166.3 


2985±267 


4.28±45 



Table 1. Emission-line properties in the integrated spectrum of SW022550 Column (1) - Line ID. Column (2) - Rest-frame wavelength 
in A. Column (3) - Redshift. Column (4) — Observed wavelength in A. Column (5) — Observed line width (Full Width at Half Maximum) 
in A, not eorrected for instrumental resolution. Column (6) — FWHM corrrected for instrumental resolution in km s~^. Integrated 
emission-line flux in 10~^® erg cm~^. 



Line ID 


Ao 


z 


Aobs 


FWHM 


FWHM 


flux 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


[OIII] 


5007 


3.42312±0.0009 


221466±4 


95.9±6 


1275±77 


40.9±2.5 


[OIII] 


4959 


3.44239±0.0005 


219380±5 


93.9±6 


1259±82 


14.6±1.0 


H/3 


4861 


3.42438±0.0009 


215069±5 


69.1±5 


929.6±70 


6.8±0.6 


[Oil] 


3727 


3.42480±0.0009 


164912±3 


105.1±7 


1895.1±120 


19.4±0.1 



Table 2. Emission-line properties in the integrated spectrum of SW022513 Column (1) - Line ID. Column (2) - Rest-frame wavelength 
in A. Column (3) - Redshift. Column (4) — Observed wavelength in A. Column (5) - Observed line width (Full Width at Half Maximum) 
ill A, not corrected for instrumental resolution. [Oil] was fitted with a single component at a rest-frame wavelength of 3727A. Column 
(6) — FWHM corrrected for instrumental resolution in km s~^. Integrated emission-line flux in 10~^^ erg s~^ cm~^. 



Line ID 


Ao 


z 


Aobs 


FWHM 


FWHM 


flux 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


[OIII] 


5007 


3.4227±0.0009 


22144±4 


46±3 


579±40 


1.22±0.09 


[OIII] 


4959 


3.4227±0.0010 


21932±5 


45±5 


579±70 


0.40±0.06 


[OIII] 


5007 


3.4032±0.0018 


22047±9 


374±32 


5087±447 


2.58±0.25 


[OIII] 


4959 


3.4031±0.0058 


21835±28 


370±95 


5087±1309 


0.85±0.25 


H/3 


4861 


3.4247±0.0010 


21504±5 


26±5 


369±40 


0.32±0.05 


H/3 


4861 


3.422 ±0.0010 


21491±5 


72±15 


1000±73 


0.30±0.07 


[Oil] 


3727 


3.4240±0.0010 


16488±3 


58±4 


1037±73 


1.19±0.09 


[Nelll] 


3869 


3.4248±0.0010 


17117±4 


41±4 


682±82 


0.37±0.06 



Table 3. Emission-line properties near the continuum peak of SW022513 extracted from a 0.4" X 0.4" aperture, including the broad and 
narrow components of [OIII]. Column (1) - Line ID. Column (2) - Rest-frame wavelength in A. Column (3) - Redshift. Column (4) 
— Observed wavelength in A. Column (5) — Observed line width (Full Width at Half Maxinmm) in A, not eorrected for instrumental 
resolution. Column (6) — FWHM corrrected for instrumental resolution in km s~^. Integrated emission-line flux in 10~^^ erg cm~^. 
We adopt the redshift of the narrow H^ component as systemic redshift for SW022513 (see text). 



Line ID 


Ao 


z 


Ao6s 


FWHM 


FWHM 


flux 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


[OIII] 


5007 


3.4231±0.0009 


22146ib4 


37±2 


455±29 


3.34±0.22 


[OIII] 


4959 


3.4231±0.0009 


21934±5 


37±3 


455±38 


l.lOibO.l 


[OIII] 


5007 


3.4143±0.0016 


22102±8 


159±19 


2152±257 


1.34±0.2 


[OIII] 


4959 


3.4143±0.0059 


21890±29 


158±66 


2152±907 


0.44±0.23. 


H/3 


4861 


3.4238±0.0010 


21505±5 


43±5 


553±70 


0.64±0.1 


[OH] 


3727 


3.4245±0.0009 


16490±4 


42±3 


734±65 


0.93±0.1 


[Nelll] 


3869 


3.4254±0.0012 


17121±5 


34±8 


559±136 


0.27±0.08 



Table 4. Emission-line properties near the peak of the [OIII]A5007 line emission in SW022513 extracted from a 0.4"x0.4" aperture, 
including the broad and narrow components of [OIII]. Column (1) - Line ID. Column (2) - Rest-frame wavelength in A. Column (3) - 
Redshift. Column (4) - Observed wavelength in A. Column (5) - Observed line width (Full Width at Half Maximum) in A, not corrected 
for instrumental resolution. Column (6) — FWHM corrrected for instrumental resolution in km s~^. Integrated emission-line flux in 10~^® 
erg cm~^. 
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Line ID 


Ao 


z 




FWHM 


FWHM 


flux 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


[OIII] 


5007 


3.4253±0.0009 


22157±5 


47±3 


604±42 


1.43±0.11 


H/3 


4861 


3.4254±0.0009 


21513±4 


20±2 


165±23 


0.29±0.05 


[Oil] 


3727 


3.4262±0.0010 


16496±4 


35±3 


606±65 


1.12±0.15 


[Nelll] 


3869 


3.4273±0.0010 


17128±4 


27±5 


409±76 


0.53±0.13 



Table 5. Emission-line properties in the elongated northern part of SW022513 extracted from a 0.7" x 0.7" aperture. Column (1) - 
Line ID. Column (2) - Rest-frame wavelength in A. Column (3) - Redshift. Column (4) - Observed wavelength in A. Column (5) - 
Observed line width (Pull Width at Half Maximum) in A, not corrected for instrumental resolution. Column (6) - FWHM corrrected for 
instrumental resolution in km . Integrated emission-line flux in 10~^® erg cm~^. 



